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ABSTRACT 

Context. Lithium is a fragile element, which is easily destroyed in the stellar interior. The existence of lithium-rich giants still repre- 
sents a challenge for stellar evolution models. 

Aims. We have collected a large database of high-resolution stellar spectra of 824 candidate thick-disk giants having 2 MASS pho- 
tometry and proper motions measured by the Southern Proper-Motion Program (SPM). In order to investigate the nature of Li-rich 
giants, we searched this database for giants presenting a strong Li I resonance line. 

/Wefhods. We performed a chemical abundance analysis on the selected stars with the MOOG code along with proper ATLAS-9 model 
atmospheres. The iron content and atmospheric parameters were fixed by using the equivalent width of a sample of Fe lines. We also 
derive abundances for C, N, and O and measure or derive lower limits on the '"C/'^^C isotopic ratios, which is a sensible diagnostic of 
the stars evolutionary status. 

Results. We detected five stars with a lithium abundance higher than 1.5, i.e. Li-rich according to the current definition. One of them 
(SPM-313132) has A(Li)>3.3 and, because of this, belongs to the group of the rare super Li-rich giants. Its kinematics makes it a 
likely thin-disk member and its atmospheric parameters are compatible with it being a 4 Mq star either on the red giant branch (RGB) 
or the early asymptotic giant branch. This object is the first super Li-rich giant detected at this phase. The other four are likely low- 
mass thick-disk stars evolved past the RGB luminosity bump, as determined from their metallicities and atmospheric parameters. The 
most evolved of them lies close to the RGB-tip. It has A(Li)>2.7 and a low ''C/'^'C isotopic ratio, close to the cool bottom processing 
predictions. 

Key words. Stars: abundances; Stars: atmospheres; Stars: chemically peculiar; Galaxy: disk 



1. Introduction 

Lithium is a fragile element, which is easily destroyed in the 
stellar interiors at temperatures higher than ~2.5xl0^ K. As soon 
as a star evolves to the red giant phase, the convective envelope 
deepens and bring to the stellar surface material that has been 
exposed to high temperatures in the stellar interior, causing an 
overall dilution of Li in a giant's atmosphere. Red giant stars are, 
therefore, expected to present a low amount of Li. Starting froni 
the meteoritic abundance of A(LiQ=3.3 dex ( Grevesse & Sauval 



|1998| l, standard models predict the dilution of Li in the giants' 
atmosphere down to a level of about 1 .5 dex, depending on the 
stellar mass and metallicity ( jIben[|T%7a|b| l. Giants having A(Li) 
exceeding this value are usually termed Li-rich and represent a 
challenge for standard stellar evolution models. 

About 1-2% of K-giants are Li -rich, according to the above 
definition ( |Brown et aL] [19891 [Smith et al.| [1995). In order to 
explain their existence, it was suggested that dilution might not 
have been effective in these stars or that the atmospheric Li abun- 
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A(Li)=log 2^ + 12 



dance might have been increased by the ingestion of a planet or a 
brown dwarf because of the expansion of the stellar atmosphere 
(Siess & Livio 1999 1. Alternatively, a Li production phase fol- 
lowing the first dredge-up dilution has been proposed. 

The possibility that Li-rich giants avoided dilution and sim- 
ply preserved their initial lithium abundance is contradicted, 
however, by the low ^^C/'^C carbon isotopic ratios measured ( da 
!Silva et al. 1995 1, which suggests that mixing has taken place in 
those stars. Moreover the existence of a few stars with Li abun- 
dances exceeding the meteoritic value ("super" Li-rich giants, 
e.g. jBalachandran et aL] |2000j ) argue against a preservation of 
the primordial Li abundance. At the same time, the ingestion of 
a planet or a brown dwarf would also lead to an increas e of the 
''Be abundance, which is not observed in Li-rich giants (jMeloet 
[all [2005] ). 

On the other hand, there is agreement among the investiga- 
tors that Li can be enhanced in the giant's atmosphere through 
the Cam eron & Fowler| ( |l97l| hereafter CF71) ^Be-transport 
mechanism. In order for this mechanism to be effective, ^He 
should first be circulated from the stellar atmosphere down to 
regions at temperatures high enough for ^He burning through 
the ^He{a, 'x)^Be reaction. Then, the produced ^Be should be 
circulated up to the stellar surface where it can decay into ^Li 
by electron capture: ^Be(e", v)^Li. In luminous (M/,o;=-6 to -7) 
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intermediate mass asymptotic giant branch (AGB) stars, the H- 
buming shell is in contact with the stellar convective envelope 
and, therefore, ^He burning takes place partially under convec- 
tive conditions (hot-bottom burning). In low-mass stars, the con- 
vective envelope and the H-buming shell are, however, not in 
touch and the existence of a mixing mechanism has to be postu- 
lated to connect material from those two zones. 

Besides the standard mixing caused by the first dredge-up, 
an additional extra-mixing episode is known to take place at the 
luminosity bump along the red giant branch (RGB) in low-mass 
stars ( |Gratton et aLl|2U04l i. [Charboniiel & Balachandran| ( |2000l 
hereafter CBOO) took advantage of Hipparcos parallaxes to ac- 
curately place Li-rich giants in the Hertzspung-Russell diagram, 
and revealed that some stars were in the process of complet- 
ing the dilution and were, hence, eiToneously classified as Li- 
rich. Besides those stars, truly Li-rich giants tend to cluster into 
two groups: at the RGB-bump luminosity in low-mass stars and 
at the early AGB in intermediate mass stars. At these phases 
the mean molecular weight discontinuity left over by the first 
dredge-up is erased, allowing the occurrence of extra-mixing 
processes. Most of the Li-rich giants present, however, '^C/'^'C 
compatible with first dredge-up predictions (CBOO). The Li-rich 
stage would then be produced by a relatively shallow extra- 
mixing episode, a precursor of the extra-mixing phenomenon 
that further lowers the '^C/'^C isotopic ratio in more evolved 
stars (see Charbo nnel & Do Nascimento|[T998l CBOO). After the 
stellar atmosphere is enriched in lithium, convection begins to 
again dilute it by exposing surface material to the higher temper- 
atures of the stellar interiors. Therefore, the Li-rich phase would 
be by nature a short-lasting phase, which agrees with the low 
number of detected stars. '^C/'-'C isotopic ratios lower than first 
dredge-up prediction may indicate that the dilution process has 
started again and Li should be lower than the peak value. 

[Charb onnel & Balachandran| (|2000|l certainly provide a 
valuable framework for the investigation of Li-rich giants. 
Nevertheless, the authors do not embrace the whole complex 
phenomenology of Li-rich giants. For instance, HD 77361 is 
a super Li-rich s tar lying clo se to the RGB -bump, but with 



Li-rich giants may rise from the standard 1-2% up to ~50% when 
rapid rotators are considered (v sin !>8 km s '). 



C/'^C=4 (Kumar & Reddy 2009), i.e. significantly lower 



than expected from first dredge-up models. Li-rich low-mass gi- 
ants were also detected along the upper RGB/ AGB in globu- 
lar clu sters (|Kraft et al.[ [l999; Smith et al.| [T999ll, the Galactic 
Bulge (lUttenthaler et al. 7007 ; Gonz alez etal.| '2009) and dwarf 
spheroidal galaxies (dSph, Dom mguez et al. 2004 ; Monaco &| 
Bonifacio] |2008 1. In particular, Dommg uez et al.| ( 2004 1 and 
Uttenthaler et al. ( 2007 1 studied a sample of C stars in the Draco 



dSph and in the Galactic Bulge, respectively. In these cases, the 
high Li abundance detected in a few stars was attributed to en- 
richment during third dredge-up episodes in thermally pulsat- 
ing AGB stars. Notice also that Li-rich stars detected in glob- 
ular clusters and in the the tidal stream of the Sagittarius dSph 
( [Monaco & Bonifacio | 2008 ) have low mass and lie either close 
to the tip of the RGB or on the AGB. As such, they do not belong 
to the second category of Li-rich giants identified by CBOO. 

Rotation has often been claimed to be the culprit for extra- 
mixing processes that happen along the RGB ( [Sweigart &[ 
Mengel 1979 1. Fekel & Balachandran ( 1993| l found that many 
rapidly rotating (v sin; =6-46 km s~'), chromospherically active 
giants present high Li abundances. They argued that the high 
observed rotational velocities are the results of momentum trans- 
ferred to the surface from the stellar interior. This transfer may 
be accompanied by the dredge-up of ^Be, whose decay increases 
the surface Li abundance. Indeed, while K-giants are usually 
slow rotators, [Drake et aL] ( [2002| ) concluded that the fraction of 



Using IRAS colors, [de La Reza et al.[ ( [19961 [1997) 1 also 
noted that many Li-rich stars presented far-infrared excess and 
speculated that the Li-rich phase was associated with a mass- 
loss episode. A few Li-iich giants certainly present evidence 
of mass loss and chromospheric activity (Balachandran et al.j 
2000t|Drake et al.[|2002| i. rFekel & Watson[ ( [l998) and ,Jasniewicz[ 



et al.[ ( 1999 1, however, analyzed samples of giants with far- 



infrared color excess and did not detect any additional Li-rich 
star. Therefore, mass loss may or may not be associated with the 
Li-rich phase. 

Several models were proposed to interpret the Li-rich gi- 
ant phenomenon. While the classical Cameron-Fowler mecha- 
nism was envisaged to take place in AGB stars (CF71) during 
He-shell flashes, the "cool bottom processing" takes place in 
RGB stars and can go on for an extended period of time as the 
star cUmbs up the RGB ( [Sackmann & Boothroyd[[1999| l. In or- 
der to be able to produce 'Li, it requires deep mixing with a 
high-speed Mp, typically Mp ^ 10^'' M^yr^^ , but often as fast 



as 1 0~'*M(7i"vr~^ A ccording to the computations of Sackmann & 
Boothroyd ( 1999 1, the ^Li in the envelope increases as time on 
the RGB passes, reaches a plateau, and then decreases again, as 
the freshly formed Li is mixed back into high-temperature re- 
gions where it is destroyed. Thus the Li-rich phase is limited in 
time. The level of Li-enrichment depends on the mixing speed 
as well as on the mixing geometry, i.e. the relative surfaces of 
upward and downward streams. 

According to'Denissenkov & Weiss ([2000') and'Denis senkov[ 
|& Herwig ( 2004) , enhanced extra-mixing may be triggered by 
the episodic ingestion of a massive planet or a brown dwarf (see 
|Carlberg et al.[[20T0) l. This could happen at any time along the 
RGB. The ingested body would naturally contribute its momen- 
tum to the stellar structure and, after the enhanced extra-mixing 
gave rise to the Li-rich phenomenon, lithium would start to be 
destroyed again by convection. This model would explain the 
existence of Li -rich giants all along the RGB at luminosities dif- 
ferent from the RGB-bump and the connection with stars with 
high rotational velocities, but would neither explain the cluster- 
ing of giants at the RGB-bump in low-mass stars nor the exis- 
tence of Li-rich giants with low rotational velocities. 

Other possible sources for the extra-mixing episodes re- 
quired to produce Li-rich stars were identified in thermohaline 
mixing ([Charbonnel & Zahn 20 07| and magneto-th ermohaline 
mixing ( [Denissenkov et al.j,2009j l. GuandaUni etal. ( [2009 here- 
after Gt'BU09) investigated the latter case and showed that 
both Li-rich and Li-poor RGB/AGB stars can be accounted 
for by models of extra-mixing induced by magnetic buoyancy. 
Their investigation does not try to interpret super Li-rich giants 
(A(Li)>2.5 dex) and still depends on several ad-hoc assump- 
tions, yet it provides a self-consistent interpretation for the ap- 
parent clustering of Li-rich stars at the RGB-bump, as well as 
the presence of Li-poor stars and the few Li-rich stars in the late 
RGB. The clump of Li -rich stars detected on the early AGB by 
CBOO were interpreted by GPBU09 as stars that actually are on 
the RGB. The first detection of a variable magnetic field in the 



Li-rich giant HD 232862 has been recently made ( [Lebre et al. 
g009l ). 

Given the complexity of the problem, it is mandatory to col- 
lect more data to better clarify the observational framework. 
Unfortunately, the intrinsic rarity of Li-rich giants requires a 
dedicated search for these kind of objects, which demands ex- 
tensive observational campaigns. 
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In the context of a survey of the Galactic thick disk, we 
constructed a large database of high-resolution, high-quality 
echelle spectra of giant s tars ([Carraro et aL 2005 Moni Bidin 



[2009 ; M oni Bidin et al.| |20l6 1. By screening this database we 



serendipitously detected the five Li-rich giants we present in the 
present work. The paper is organized as follows: in ^we outline 
the target selection criterion and present the observations and 
data reduction performed with reference to the detected Li-rich 
stars. In section ^we present the performed abundance analy- 
sis on the selected stars and in ^we present the measured target 
rotational velocities. In section ^ we use the derived chemi- 
cal abundances and atmospheric parameters together with their 
radial velocities and proper motions to investigate whether the 
kinematics of our stars are compatible with membership to the 
thick-disk. Finally, in sections ^and ^we discuss and sum- 
marize our findings. 



2. Observations and data reduction 



Our investigation is based on the sample defined by Girard et al. 
2006| hereafter G06). It consists of 1 196 red giant stars within 



5 degrees from the south Galactic Pole, with infrared 2 MASS 
photometry (Skrutskie et al. 20 06]) an d absolute proper motions 
from the SPM3 catalog ( ,Girard eTaT] [20041 . 

The sample was defined from an IR color-magnitude dia- 
gram in the color range 0.7<J-K <L1, to select intermediate 
metallicity stars, i.e. thick-disk members. Contamination by halo 
stars should account for ~8% of the selected objects (G06). The 
sample is confined to a cone perpendicular to the Galactic disk, 
which G06 estimated to be volume-complete up to at least 3 kpc 
from the Galactic plane. A sloped-cutoff limit at faint magni- 
tudes avoided contamination by dwarf stars, except for the clos- 
est ones (d<63 pc). 

High-resolution, high-quality echelle spectra were acquired 
with four different instruments for a total of 824 stars belong- 
ing to the G06 sample during 38 observing nights between 2005 
and 2007. Telescopes with dififerent collecting powers were em- 
ployed to efficiently secure spectra with roughly the same resolu- 
tion and quality for targets spanning a wide range of magnitudes 
(5<V<16). 

We inspected by eye all spectra in the database in the region 
around the lithium resonance doublet at 6707.8 A. The line was 
detected in several stars, but only the five objects in Table[2]pre- 
sented a strong line. Indeed, the Li subordinate line at 6103.6 A 
was detected in these stars only. A detailed analysis of all objects 
with detected Li lines will be presented elsewhere. We focus here 
on the five stars that alone present strong lithium lines. 

Stars discussed here were observed either with the 
Echelle/duPont or the FEROS/2.2m spectrographs. Below we 
briefly summarize the data reduction performed with reference 
to these observations. Details can be found in |Moni Bidin| ( |2009l l 
and Mo ni Bidin et al.| ( |2010| l. In Table[T] we present the SPM 
identification number (which will be used henceforth), 2 MASS 
designations, date of the observations, instruments, observatory, 
exposure times, and resolutions for our stars. Table|2]reports the 
target coordinates, the 2 MASS infrared magnitudes, the SPM3 
proper motions and the reddening as derived from the [Schlegel 
|et^([l998jl maps. 

During each run we observed several bright standard stars 
with accurate parameters (radial and rotational velocities, fun- 
damental parameters and metallicity) available in the literature. 
Among them, the well known Li-rich giant HD 787 ( jCastilho et 
[ak] |2000l was also observed. This star is analyzed here as well. 



Data reduction was perfoimed using standard IRAlj^ tasks 
and included bias, flat-fielding, order tracing, background and 
sky subtraction, and extraction. Standard Th-Ar arcs were used 
for the wavelength calibrations. For the duPont spectra we 
used "milky" flats, obtained with a diffusion filter set be- 
fore the CCD. These frames provided the needed wavelength- 
independent pixel-to-pixel response variation but, at odds with 
standard spectroscopic flat-field frames, the "milky" flats did not 
correct the fringing pattern that affected our spectra longward of 
about 7500 A. For these, we created a fringing image from the 
Dome Flat, reduced as a normal spectrum and then normalized, 
and we divided all science spectra by this image in the reddest 
orders, where the fringing pattern became visible. The resulting 
correction was deemed satisfactory. 

The radial velocities (RVs) of the observed stars were mea- 
sured with a cross-correlation (CC) technique ( Tonry & Davis| 
[1979) by means of the fxcor IRAF task. Radial velocity stan- 
dards acquired with the same instrument were used as templates 
for the CC. After the radial velocity measurements, all spec- 
tra were shifted to the laboratory wavelength. The IRAF task 
rvcorrect was used to calculate the Earth motion and convert 
the observed RVs to the heliocentric system. Measured RVs are 
reported in Table[3]as well as the spectral signal-to-noise (S/N). 

3. Chemical abundance analysis 

For each target we obtained an initial estimate of the stellar effec- 
tive temperature (Tss) from the dereddened infrared (J-K) color 
along with the Alonso et al. ( 19991) calibrations. The reddening 



was derived from the COBE reddening maps (Schlegel et al. 
1998), obtaining E(7 - K) from E(B - V) by means of the trans- 
formations of [Cardelli et al.|( 1989). The surface gravity (logg) 



was derived with the aid of theoretical isochrones ( Marigo et al. 



'2008), adopting an age of lOGyr (see also IFeltzing & Bensby 
2009), a metallicity Z=0.004, and a distance estimated from a 
photometric parallax, by assuming that target stars followed the 
same (J-K) vs K mean ridge line as the globular cluster 47 Tuc. 
Indeed, the thick disk age and metallicity are similar to those of 



this cluster ( Wyse & Gilmore 2005 i. The microturbulent veloc- 
ity ^ was evaluated from the derived gravity adopting the scahng 



2008 



relation: ^ (km s ')= -0.254 x logg H- 1.930 ( [Marino et al. 
hereafter M08). 

These initial parameters were adopted to compute a proper 
model atmosphere for each star using the the Linux port of the 
A TLAS -9 code ([Kurucz) [1993a , Sbor done, Bonifacio, CasteliT] 
|& Kurucz||2004[ ). Equivalent widths (EWs) of a selected sample 
of iron lines (see M08) were measured with a Gaussian fitting. 
For each star, the measured EW s and the model atmosphere were 
then used within the M00G(^ ( [Sneden[ [l973| local thermody- 
namic equilibrium (LTE) stellar line analysis code to refine the 
atmospheric parameters. In particular, the stellar Tetj was fixed 
by imposing that the abundance derived for each line should 
be independent of the transition excitation potential. ^ was es- 
timated by minimizing the dependence of the abundances from 
the measured EWs, and log g by requiring that Fe I and Fe II lines 
gave the same iron abundance, within the errors. The stellar at- 
mospheric parameters eventually adopted are reported in Table[3] 
together with the derived iron abundance. 



- IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 1. Log of the observations for the program stars. The SPM and 2 MASS identification numbers are also reported together 
with the spectral resolution of the obtained observations. 



SPM 


2MASS 


Date 


Exposure Time 
(s) 


A|^A 


Instrument 


Telescope 


Observatory 


142173 


00321256-3834022 


2006-09-13 


1800 


38 000 


Echelle 


2.5m Irenee du Pont 


Las Campanas 


171877 


00392022-3131354 


2006-08-23 


900 


48 000 


FEROS 


2.2m ESO-MPI 


La Silla 


225245 


00544638-2735304 


2006-09-11 


2400 


38 000 


Echelle 


2.5m Irenee du Pont 


Las Campanas 


313132" 


01202066-3409541 


2006-09-12 


300 


38 000 


Echelle 


2.5m Irenee du Pont 


Las Campanas 


343555* 


01294200-3015464 


2006-09-12 


850 


38 000 


Echelle 


2.5m Irenee du Pont 


Las Campanas 


HD 787^ 


00120998-1756177 


2006-08-20 


100 


48 000 


FEROS 


2.2m ESO-MPI 


La Silla 



" Alternative ID: CD-34 510. 
* High-velocity star. 
^ Alternative ID: HR37. 



Table 2. Coordinates, photometry, reddening and proper motions (pm) of the program stars. 



SPM 


a(J2000) 


(5(J2000) 


1 


b 


K 


J-K 


E(B-V) 


E(J-K) 


pm(RA) 


pm(Dec) 








deg 


deg 










mas/yr 


mas/yr 


142173 


00:32:12.56 


-38:34:02.3 


321.099 


-77.874 


9.29 


0.78 


0.018 


0.009 


-1.2±2.0 


-19.4±1.6 


171877 


00:39:20.23 


-31:31:35.5 


333.144 


-84.872 


8.12 


0.91 


0.019 


0.010 


6.8±1.6 


-7.7±2.5 


225245 


00:54:46.38 


-27:35:30.4 


245.156 


-89.126 


8.79 


0.88 


0.014 


0.007 


4.5±1.2 


-3.3±1.1 


313132 


01:20:20.66 


-34:09:54.1 


263.253 


-80.616 


7.04 


0.73 


0.031 


0.016 


26.4±3.1 


-5.2±2.4 


343555 


01:29:42.01 


-30:15:46.4 


235.664 


-81.047 


8.46 


0.72 


0.017 


0.009 


27.9±1.7 


-25.3±1.3 


HD787 


00:12:09.99 


-17:56:17.8 


76.321 


-77.101 


1.85 


0.85 


0.029 


0.015 







Table 3. Program stars atmospheric parameters. The spectra signal-to-noise ratios are also indicated, as well as the measured radial 
velocities. 



SPM 




logg 




A(Fe) 


[Fe/H] 


S/N 


Vhelio 




K 




km s ' 






@670nm 


km s"' 


Sun 


5777 


4.44 


0.80 


7.50 








142173 


4330 


1.50 


1.35 


6.81 


-0.69 


118 


38.0 ±0.7 


171877 


3930 


1.10 


1.40 


6.73 


-0.77 


109 


-15.8 ±0.5 


225245 


3920 


0.65 


1.50 


6.51 


-0.99 


110 


3.9 ±0.6 


313132 


4530 


2.00 


1.20 


7.51 


-hO.01 


42 


7.7 ±0.6 


343555 


4530 


2.25 


1.00 


6.88 


-0.62 


106 


119.8 ±0.6 


HD787'' 


4000 


1.50 


1.50 


7.50 


+0.00 


>150 


-6.0 ± 0.9 


HD787 


3870 


0.85 


1.27 


7.54 


+0.04 


>150 


-6.0 ± 0.9 



" Atmospheric parameters are assumed from |Castilho et aL] ( |2000| ). 



The abundance of Li, C, N, O, and the '^C/'^C isotopic ra- 
tio were obtained with spectral synthesis of selected regions us- 
ing the MOOG code along with an ATLAS-9 model atmosphere 
computed for each star with the final adopted atmospheric pa- 
rameters. The atomic and molecular linelists were taken from 
the Kurucz databas^ 

For the hottest stars in the sample, i.e. #142173, #313132, 
and #343555, it was possible to measure the carbon abundance 
from the Ci band at 5635.5 A. For stars #171877, #225245 arid 
HD787 carbon was measured from the CH G-band at 43 10 A. 
For #343555 we managed to measure carbon from both bands 
and obtained an abundance difference of 0.05 dex, which proves 
the consistency of the two measurements. Oxygen was measured 
by synthesis of the fine at 6300 A. The Nil fine at 6300.34 A 



http : //kurucz . harvard . edu/ 



was taken into account in the synthesis of the Oxygen line. For 
the Ni I line we adopted the log gf value derived by | Johansson] 
et al. ( 2003| l, namely loggf=-2.11. The nitrogen abundance was 
derived by synthesis of the 8001-8005 A spectral region, which 
is rich in CN lines. The same region was also used to set the 
'^C/'^C isotopic ratio, thanks to presence of the '^CN line at 
-8004.7 A and the '^CN lines at 8003-8004 A. Owing to the low 
residual intensity of the '^CN line, for stars #142173, #171877, 
and #343555 we provide only lower limits to the carbon isotopic 
ratio. 

Because the C-N-O elements are partly bound together in 
molecules for stars cooler than ~4500 K, their abundances were 
derived together in an iterative process. Solar abundances for Fe, 
C, N, and O were obtained in the same way using the Kurucz 

[M4| l. The C2 



solar spectrum for the purpose (Kurucz et al 



band at 5635.5 A is not visible in the solar spectrum, so we used 
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Fig. 1. Sample of the target stars' spectra in the region of the Li I Fig. 2. Sample of the target stars' spectra in the region of the Li I 



resonance doublet. 



subordinate line. 



the CH G-band. The derived solar abundances are reported in 
Table m 

Afterward, we measured the lithium abundances from both 
the resonance doublet at 6707.8 A and th e subordinate line at 
6103.6 A 



For the former, we adopted the Reddy et al.| ( |2002 i 
hyperfine structure. We also added to the linelist the two Fe I and 
Til hnes at 6707.433 A and 6708.125 A , respectively, listed in 
Gonzalez et al. ( |2009| hereafter G09). These lines are not present 
in the Kurucz database^ but are available at the Vienna Atomic 

5 



Line Database (VALEjj |Kupka et al. 2000) ) and are required to 
obtain a good fit of the Li I doublet at 6707.8 A. For these lines 
we adopted the oscillator strength (log gf) calibrated on Arcturus 
by G09. 

Figures [T] [2] and |3] present the spectra of our stars and 
HD 787 around the two Li I Unes and in the 8001-8005 A region. 
For each star, the best-fitting synthetic spectrum (continuous 
line) is superimposed on the observed one (filled squares). The 
measured C, N, O, and Li abundances are reported in Table?] 
For star #343555, the lithium abundance was measured from the 
6707.8 A line, while we derive an upper limit only from the sub- 
ordinate line, which is barely detected in the observed spectrum. 

The first of the two HD 787 entries in Table[4] reports the Li 
abundances we obtained adopting the atmospheric parameters in 
Castilho etal] ( [20001 see first HD 787 entry in Table[3]l. We de- 
rived Li abundances of 2.03 and 2.14dex for the 6707.8 A and 



Castilho et 



6103.6 A lines, respectively. This agrees well with^ 

al. ( 2000| l, who measured 2.0 and 2.2 dex (LTE analysis). These 
authors did not measure the C, N, and O abundances and, there- 
fore, we adopted the Melendez et al. (2008) values, which are 
also reported in the Table. We also report the ^^C/'^C isotopic 
ratio measured by |da Silva et al.| ( [l995| l. 

For HD 787, however, we derive a slightly cooler temper- 
ature (see second HD787 entry in Table[3]) with respect to 
|CastiUio et al.| ( [2000) i. We therefore performed an independent 




8002 B003 8004 8005 8002 8003 8004 8005 

Wavelength (A) 

Fig. 3. Sample of the target stars' spectra in the region used to 
measure the N abundance and the '^C/'-'C isotopic ratio. 



abundance analysis, following the same procedure as outlined 
above. The second HD 787 entry in Table|4] reports the C, N, O, 
and Li abundances we derived for HD 787 with our analysis. The 
'^C/'^C isotopic ratio we measured also perfectly agrees with|da| 
ISilva et al.|(|l^. 



Lind et al. ( 2009 | l calculated the coiTection to the Li abun- 



http : //vald . astro . univie . ac . at/~vald/php/vald . php 



dances measured from both the resonance and the subordinate 
lines at 6707.8 A and 6103.6 A due to non-LTE eff^ects (NLTE). 
We used their grid to evaluate the correction to apply to our 
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Table 4. Measured abundances for the program stars. The measured stellar rotational velocities are also indicated. 



SPM 


[Fe/H] 


[C/Fe] 


[N/Fe] 


[O/Fe] 


C/O 


'-c/"c 


A(Li) 


A(Li) 




A(Li);vLr£ 


vsin; 


















670.8nm 


610.3nm 


670.8nm 


610.3nm 


km s" 


1 


Sun 


7.50 


8.49 


7.95 


8.83 


















142173 


-0.69 


+0.02 


+0.02 


+0.20 


+0.30 


>15 


2.80 


2.52 


2.73 


2.72 


2.2 ± 


0.7 


171877 


-0.77 


-0.23 


+0.42 


+0.34 


+0.12 


>15 


2.49 


2.29 


2.46 


2.52 


<1.8 




225245 


-0.99 


-0.44 


+0.54 


+0.08 


+0.14 


8 


2.90 


2.71 


2.83 


2.93 


<1.1 




313132 


+0.01 


+0.07 


+0.01 


-0.07 


+0.63 


13 


3.45 


3.35 


3.37 


3.56 


3.3 ± 


0.5 


343555 


-0.62 


+0.09 


-0.08 


+0.19 


+0.36 


>10 


1.79 


<1.94 


1.94 


<2.11 


2.5 ± 


0.8 


HD 787" 


+0.00 


-0.24 


+0.33 


+0.02 


0.21 


15 


2.03 


2.14 


2.09 


2.40 






HD787 


+0.04 


-0.38 


+0.09 


-0.22 


0.27 


15 


1.77 


2.12 


1.98 


2.40 







" NLTE corrections are calculated according to Lind et al. 
for #171877, #225245, and HD787 (2nd instance); logg= 
instance). 



1 2009 1. Owing to the grid boundaries, we assume: ^=2.0km s^'for all stars; reft=4000K 
:1.0 for #225245 and HD787 (2nd instance); [Fe/H]=0.0 for #313132 and HD787 (2nd 



' Atmospheric parameters are assumed from Castilho et al. 1 2000 1; C, N and O abundances are taken from Melendez et al. 1 2008 i. Carbon isotopic 



ratio is taken from ,da Silva et al.,^1995^ . 



measures. In columns 10 and 11 of Table[4] we report the Li 
abundances coiTected for NLTE. These values have to be re- 
garded only as referential. Indeed, owing to the grid boundaries, 
we had to adopt in the computation ^=2. 0km s^'for all stars; 
reff=4000K for #17 1877, #225245, and HD787 (2nd instance); 
logg^l.O for #225245 and HD787 (2nd instance); [Fe/H]=0.0 
for #313132 and HD 787 (2nd instance). 

We finally estimated the eiTors on absolute abundances and 
abundance ratios owing to the noise of the spectra and to the 
uncertainties on atmospheric parameters. For the '^C/'^C ratio 
only the error owing to the spectral S/N is significant because a 
change in temperature, gravity, or microturbulence does not al- 
ter the relative intensity of the '^CN or '^CN molecular lines. For 
Fe, the uncertainty owing to the noise is simply the eiTor on the 
mean as given by MOOG. For the other elements, whose abun- 
dances were obtained by spectrosynthesis, the errors were eval- 
uated by comparing observed spectra with synthetic ones calcu- 
lated with a step of 0.05 dex in abundance or 2 in i^C/i^C. We 
estimated by eye the abundance interval that is still compatible, 
within the noise, with the observed spectrum. The error owing to 
the noise was assumed to be the semi-amplitude of this interval. 

As for the atmospheric parameters, uncertainties were esti- 
mated by changing one parameter at the time and measuring the 
abundances again. The error for a given parameter is given by 
the difference between the old and the new abundances. We as- 
sumed ±100 K, +0.2 dex, and ±0.10 km/s as the typical eiTors 
for Teff, logg and ^, respectively, and analyzed stars #313132 
and #225245, which cover the entire Teff range of our sample. 
Results are reported in Table[5] The last column in the table re- 
ports the total error, which is defined as 



0"tot 



4 



+ erf , , + cri. 

log(g) f 



4. Rotational velocities 

The projected rotational velocity of the target stars, vsin/, was 
measured from the width of the auto-cross correlation function 
(CCF) of their spectrum: each one was cross-correlated (Tonry 
|& Davis|[T979if with itself through the fxcor IRAF task, and the 
peak of the resulting function was fitted with a Gaussian pro- 
file to derive its width cr. We verified that the resulting width 
changed when correlating different spectral ranges, but only by 
the addition of a color-dependent constant on all standard and 
target stars. This did not affect the measurements, since it only 



implied a different value for the calibrated constants of the cr- 
vsin/ relation. We finally adopted the range 6000-6500 A, a 
high-S/N region free of telluric bands and strong stellar lines, 
avoided because of their non-Gaussian wings. 

The square of the width of the auto-CCF (cTobs) can be ex- 
pressed as the quadratic sum of the rotational broadening (cTj-oj) 
and a term including all the causes of line broadening other than 
rotation (ctq). The relation between the measured (Tabs and the 
projected rotational velocity can thus be written as 



o-^„t ^k-(v sin if = o-f 



obs 



(1) 



where kisa constant coupling v sin / to the line rotational broad- 
ening, ctq is the sum of an instrumental constant and another 
term depending on the physical conditions in the stellar atmo- 
sphere (temperature, gravity, macro- and micro-turbulence). As 
a first approximation, ctq can be expressed as a function of color 
d Benz & Mayo r| |1984j [Queloz et aL] [I998)l. Previous investi- 
gations showed that a linear relation between ctq and (J - Kf 
is the most appropriate ( |Melo et"aL |2001) but, having at our 
disposal only five red giant standard observed for the calibra- 
tion of the coefficients, we preferred the lower-order expression 
(Ty - a + b ■ (J - Kf. This approximation is allowed because 
our target stars are distributed in a narrow range of colors, and 
the systematics caused by a lower-order fit are expected to be 
negligible. In summary, renaming the involved constants appro- 
priately, the relation between the measured o-gbs and v sin / can 
be written as 



(v sin if = A ■ cr^^^ +B + C(J- Kf. 



(2) 



The coefficients (A,B,C) were calibrated through a two- 
variables least-square procedure, by means of the five red gi- 
ants standard observed during the same runs as our targets, with 
(J - K)-Q.5Q-Q.2t5 and projected rotational velocity from the 
catalog of [Glebocki & Stawikowski| (|2000|. The FEROS and 
duPont runs were calibrated separately, because of their different 
instrumental broadening. 

The validity of the solution was tested comparing the v sin ; 
of the standard stars calculated through Equation [2] with its liter- 
ature value. The results are given in Table|6] The mean difference 
is negligible and the dispersion of the differences with literature 
values is 0.6 km s"'. 

The rotational velocities measured for our program stars are 
given in Table |4] The errors were estimated varying (J - K), 
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Table 5. Estimated errors on absolute abundances and abundance ratios for the hottest (#313132) and coolest (#225245) stars in our 
sample. 



A(E1) 


S/N 


ATeff=100 


Alog(g)=0.20 


Af=0.10km/s 


0"tot 


A[Fe/H]#3i3i32 


0.02 


+0.04 


+0.02 


-0.05 


0.07 


A[Fe/H]#225245 


0.01 


+0.01 


+0.03 


-0.05 


0.06 


A[C/H]#3i3i32 


0.03 


-0.02 


+0.04 


+0.00 


0.05 


A[C/H]#225245 


0.05 


-0.03 


+0.04 


+0.01 


0.07 


A[N/H]#3i3i32 


0.05 


-0.01 


+0.12 


-0.01 


0.13 


A rXT /TUT 


A AO 

O.Os 


-U.IU 


+U.i J 


+U.UU 


n 1 o 


^L^/nJ#3i3]32 


u.uj 


+U.Ui 


-lO no 


+U.UU 


n 1 n 


A r/^ /TUT 

A[0/HJ#225245 


0.02 


+U.UZ 


+U.UO 


+U.UU 


n m 
U.U/ 


Al, L,;#3|3i32 


c 

J 


n 
U 


u 


U 


c 

J 


A / 12/^ / 1 3 f~^\ 
i\y W C;#225245 


2 











2 


A(l0ge Ll67o)#313132 


n 1 n 


+U. lo 


+U.Uj 


n no 
-U.U/ 


n 1 n 
U.i / 


A(log6 Li670)#225245 


0.08 


+0.15 


+0.03 


-0.05 


0.18 


A(log6 Li6io)#313132 


0.03 


+0.10 


+0.00 


+0.00 


0.10 


A(log6 Li6,o)#225245 


0.02 


+0.13 


-0.01 


+0.00 


0.13 


A[C/Fe]#3i3i32 


0.03 


-0.06 


+0.02 


+0.05 


0.09 


A[C/Fe]#225245 


0.05 


-0.04 


+0.01 


+0.06 


0.09 


A[N/Fe]#3i3i32 


0.05 


-0.05 


+0.10 


+0.04 


0.13 


A[N/Fe]#225245 


0.08 


-0.11 


+0.10 


+0.05 


0.18 


A[0/Fe]#3i3,32 


0.05 


-0.03 


+0.07 


+0.05 


0.10 


A[0/Fe]#225245 


0.02 


+0.01 


+0.03 


+0.05 


0.06 



o"obsj and the input quantities of the calibration (color, vsin/, 
and CCF width of the standard stars), by ±lcr. The error on 
iJ - K) was taken from the 2MASS catalog, while an uncer- 
tainty of 0.05 km s~' was assumed for the measured cr, accord- 



ing to Equation 6 of Melo et al. ( 2001j ). For very slow rotators, 
Cobs ~ ctq and measurement errors can cause (Tobs < o"(), thus 
leading to an unphysical negative rotational velocity. This hap- 
pened for two target stars. The u pper limit of thei r vsin/ was 
estimated as v sin in,ax - A V2croe ( Melo et al. 2001 1, where e is 
the uncertainty on ctq, evaluated with the same procedure used 
to derive the eiTors on rotational velocity. 

The derived projected rotational velocity is very low for all 
stars. This general result is confirmed by a simple inspection of 
the CCFs, because the o-gbs of the targets never exceeded the 
analogous quantity for more than a couple of standard stars. 
Although the width also depends on the spectral type, this indi- 
cates that V sin / of all the stars must be in the range spanned by 
the standards, i.e. 0-5 km s"'. The resulting eiTors are very sim- 
ilar to the dispersion of the differences with the Uterature value 
for the standard stars. 



5. Thick-disk membersliip 

The stars under analysis are drawn from a sample of likely thick- 
disk stars, as from the target selection criterion (see Further 
insights into the actual nature of stars in our sub-sample can now 
be obtained by coupling the star kinematics to stellar isochrones 
and the estimated metallicities and atmospheric parameters. 

We estimated the star distances as follo ws. For each star, we 
first generated an appropriate isochrone|^ ( jMarigo et al. |2008| l 
using the metallicities in Table[4] An age of 10 Gyr was adopted 
for all stars at this stage ( Feltzing & Bensby{ |2009 1. Then, we 
inspected the corresponding isochrone searching for the best 
match of both gravity and effective temperature. This provides 
an estimate of the absolute magnitude which, together with the 
apparent magnitude and reddening in Table[2j allows to derive 



the star distances. Allowing for an uncertainty in the age of 
±2 Gyr results in a typical error in distance of about 10%. The 
derived distances are listed in Table|7] For further reference, no- 
tice that a variation of AreflF=±100K and Alogg=+0.5 would 
result in a variation on the derived distances of about 4% and 
17%, in the case of #142173 taken as representative of the stars 
under investigation. 

To test the consistency of the possible thick-disk member- 
ship of the target stars, we explored their kinematical and or- 
bital properties. Using the star distances in Table[7]and the star 
radial velocities (Table[3]l and proper motions (Table|2]), we es- 
timated the Cartesian Galactic coordinates and velocity compo- 
nents (see Tab le[7]i fo llowing the method described by Johnson] 
& Soderblom (1987i. Finally, we calculated the Galactic orbit 
of each star adopting a gravitational potential including the bar, 
deriving estimates for the eccentricity (e), the apo- and peri- 
galacticon (Ra and Rp), and the absolute value of the maxi- 
mum height above the plane Xmax, which are also indicated in 
Ta ble[7| Orbits we re computed following the method described 



in |Magrini et al.| ( |2010| l and [Jilkova et al~| ( |2011[ ), which the 



http://stev.oapd.inaf.it/cgi-bin/cmd 



reader is referred to for details. 

Based on the above analysis, we can draw the following con- 
clusions: 

- according to the derived orbital parameters, all stars but 
#313132 look like typical thick-disk stars; 

- #313132 has a low eccentricity and does not reach high Z. 
Besides, its orbit is confined within the solar vicinity. It 
is, therefore, reasonable to consider it as a thin-disk star. 
We note also that #313132 presents a solar metallicity (see 
Table[4]|, which is typical of thin-disk stars. 

6. Discussion 

The abundances measured for the target stars are summarized 
in Table[4] All five stars have lithium abundances higher than 
1.5 dex, and are, therefore, to be considered as Li -rich. Three 
of them present abundances exceeding 2.2 dex and one, namely 
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Table 6. Projected rotational velocity derived for the standard stars. The literatiire v sin /, and the derived values for the FEROS and 
DuPont runs, are given in columns 2-4. In the last column the Uterature references are given: 1: |de La Reza & Da Silva ( 1995 i; 2: 

' - - ■ - - ' - — - -7. 



(1982 



de Medeiros & Mayor ( 1999); 3:lde Me deiros et al.| (| T996|); 4:|Fekel & Balachandran](|T993]l; 5: |Fekel| ( |1997| l; 6: [Gray 

Pasquini et al. ^200^; 8: Smith & Dominy '{ 1979| l; 9: Gray & Tonerl 1986). When more than one estimate was available, the value 



in column 2 is their weighted average, with their standard deviation as associated error. 



ID 


V sin /(lit.) 


V sin /(FEROS) 


V sin /(duPont) 


ref. 




km s"' 


km s"' 






HDIST 


2.1±1.0 


1.3 


1.5 


1,2,3,4 




2.1±1.0 




1.6 




HD 31767 


1.1±1.3 


2.0 


2.3 


2 


HD 180540 


4.8±0.3 


4.7 


4.8 


9 


HD4128 


3.5±04 


3.6 


3.4 


5,6,7,8 


HD 15453 


1.0±1.0 


0.9 


1.3 


2 



" For HD787 only the FEROS spectrum obtained in the run reported in Table[T]was used for the abundance analysis. However, the star was 
repeatedly observed during the various duPont and La Silla runs. We report here the rotational velocities measured for the different spectra. 



Table 7. Derived distances and orbital parameters of the program stars. 



SPM 


Re 


X 


Y 


Z 


U 


V 


W 


e 










kpc 


kpc 


kpc 


kpc 


km s"' 


km s"' 


km s"' 




kpc 


kpc 


kpc 


142173 


3.6±0.3 


7.92 


0.47 


-3.48 


-195±36 


41±39 


36±9 


0.83±0.02 


3.42±0.19 


1.13±0.01 


12.41±0.14 


171877 


3.1±0.3 


8.25 


0.13 


-3.13 


17±28 


-82±37 


27±3 


0.61±0.01 


2.19±0.06 


2.00±0.01 


8.27±0.06 


225245 


7.3±0.7 


8.55 


0.10 


-7.28 


56±41 


-52±43 


6±1 


0.76±0.03 


5.11±0.22 


1.60±0.02 


11.66±0.22 


313132 


0.7±0.1 


8.51 


0.12 


-0.70 


50±11 


-165±11 


12±2 


0.28±0.01 


0.67±0.01 


5.10±0.01 


9.09±0.01 


343555 


1.1±0.1 


8.60 


0.14 


-1.07 


33± 9 


-29±20 


-84±3 


0.88±0.04 


1.93±0.10 


0.46±0.03 


7.37±0.15 




In Fig. [4] we plot our stars in the TefF vs log g plane, together 
' "~ ' ""^^ ^■2008)span- 



5000 



4000 
(K) 



3000 



Fig. 4. Target stars are plotted in the log g vs T^ff plane together 
with lOGyr theoretical isochrones of metallicity similar to the 
target stars. The symbol size is scaled according to the stars' Li 
content. 



#313132, belongs to the group of the rare super Li-rich giants 
(see, e.g., Ba lachandran et al.| [2000 [Domfnguez et al.] |2004[ 
Reddy & La mbert[|2005[ [Monaco & Bonifacio||2008[|Ktjmar & 
Reddy,|2009| i. 



with lOGyr theoretical isochrones ( Marigo et al.| 
ning the metallicity range covered by our stars. Though they 
likely cover a range of ages, Feltzing & Bensby| ( [2009 ) con- 
cluded that thick-disk stars are older than the Sun with a mean 
age of lOGyr. The position of the RGB -bump and RGB -tip are 
marked in the figure for reference. For younger ages the RGB- 
bump position on the plane moves to cooler temperatures and 
lower gravities, but for ages as young as 3 Gyr, its location lies 
still within the shaded area in the figure. 

As discussed in section ^ #313132 is likely to be a thin-disk 
star, and we will discuss its case separately. The other four stars 
seem to be more evolved than the RGB -bump, but lie below the 
RGB-tip. Their position in Fig.[4]would suggest they might have 
already suffered the extra-mixing process known to take place at 
the RGB-bump. While for stars #142173, #171877 and #343555 
we can only derive lower limits to the '^C/'^C isotopic ratio, 
star #225245 indeed presents a value lower than first dredge- 
up prediction for low-mass stars. If the Li-enrichment phase in- 
deed occurs at the RGB-bump (CBOO), this star may be in the 
process of diluting again their surface Li abundance and would 
be, therefore, similar to the HD 148293 ('2C/'^C=16, A(Li)=2) 
and HD 183492 cases C^C/^^C^9, A(Li)=2) discussed by CBOO. 
Following CBOO, if it is indeed diluting its surface Li, its posi- 
tion in Fig. |4] and its carbon isotopic ratios might provide con- 
straints on the time-scale for the extra-mixing process to occur. 
At variance with HD 148293 and HD 183492, however, #225245 
is significantly more evolved than the RGB-bump and have high 
lithium abundances (A(Li)>2.7). 

Cool-bottom processing (CBP) could still be considered as a 
viable mechanism to produce the lithium abundances observed 
in these four stars. In particular, #225245 has quite a high lithium 
abundance (A(Li)>2.7). It lies close to the RGB-tip and has quite 
a low carbon isotopic ratio ('^C/'^C=8), similar to the predic- 
tions made by [Boothroyd & Sackmann| ( |1999|| l for low-mass stars 
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^^^A(Li)=3.35; [Fe/H] = + 0.01 

Q — A(I.i)-a.52/2.7I; 

- [Fe/H|=-0.69/-0.g9 

A(l.i) = 2.29, [Fe/H]^-0,77 




5000 



4000 

Tct (K) 



3000 



Fig. 5. Target stars are plotted in the log g T^g plane together 
with theoretical evolutionary tracks of metallicity similar to the 
target stars and different masses. The symbol size is scaled ac- 
cording to the stars' Li content. 



at the RGB-tip using a model of deep mixing and the associated 
CBP. These models were proved to be effective in producing Li- 
rich stars all along the RGB ( Sackmann & Boothroyd 1999 1 up 
to very high abundances, A(Li)=4.0. Considering its low metal 
content, #225245 may also be similar to the Li-rich giants dis- 
covered in the globular clusters M 3 and NGC 362 ( ^Kraft et al.| 
[T999i [Smith etaLl[T999] l. 

Together with our stars, we plot in Fig.[5]stellar evolutionary 
2008 , 2009) of metallicity appropriate to 
cannot be consid- 



tracks (Bertelli et al. 



our stars. While the arguments in section 
ered as proof that they are thick-disk members, the locus occu- 
pied in the plane by the four stars we just discussed is consistent 
with them indeed being low-mass objects and, therefore, with 
the conclusions we have drawn above. We note that CBOO has 
shown that Li-rich giants cluster either at the RGB-bump or at 
the early AGB, in the case of low and intermediate-mass stars, 
respectively. Our stars, however, do not belong to either of these 
two categories. 

It is clear from Fig.|5]that the atmospheric parameters of star 
#313132 are perfectly compatible with a 4Mo star. Its location 
in the plane is compatible with either being on the RGB or the 
AGB. This object may accordingly belong to the second group 
of Li-rich stars identified by CBOO, i.e. stars on the early AGB, 
where the deepening of the envelope has not yet reached its max- 
imum penetration and an extra-mixing process is supposed to oc- 
cur to drive ^He from the envelope to a region hot enough for the 
Cameron & Fowler ( I97I i mechanism to take place. In this case. 



it would be similar to HD787 and HD 30834, which also have 
similar '^C/'^C isotopic ratios, compatible with first dredge-up 
predictions for similar masses (see discussion in CBOO). Star 
#313132 is, then, the first super Li-rich giant belonging to this 
group. Indeed, stars on the early AGB present on average lower 
Li abundances than stars at the RGB-bump. Guandalini et al. 
( |2009| l have interpreted this group of stars as actually lying on 




8560 6565 8560 6585 

Wavelength (A) 

Fig. 6. Sample of the target stars' spectra in the region of the H„ 
line. 



the RGB. Their models are unable to interpret super Li -rich stars 
though. We notice that it is in the 4-8 M© mass range that hot- 
bottom burning leading to Li production starts to became effec- 
tive. This is known to occur at bolometric magnitudes Mhoi - 
-7 ( [Sackmann & Boothroyd| [T992| 



Smit h et al.' 1995 

h et al. , 2008 , 2009] 

Marigo et al. 2008) suggest instead that atmospheric parameters 



-6 to 

Theoretical tracks and isochrones ( [Berte 



9} 



similar to that of star #313132 are compatible with significantly 
fainter luminosities only, namely M^oi --2.1 to -1.4, for the 
following combination of parameters: (Teff ;log g;mass)=(4400- 
4600K;1.8-2.2;3.5-4.5Mq). Therefore, the high Hthium abun- 
dance of #313132 should not have been reached through hot- 
bottom burning. 

All our stars but #343555 have Li abundances that exceed or 
are close to the maximum value allowed by the GPBU09 model 
for stars at the RGB-bump (A(Li)^2.5 dex), while more evolved 
stars would present even lower values. 

As discussed in section ^ the fraction of Li-rich giants 
seems to be highly enhanced among rapid rotators ( [Drake et"aL 
2002). The rotational velocities measured for our stars are re- 
ported in Table|4] All our targets are slow rotators, with only the 
super Li -rich star #313132 presenting a slightly higher value, 
namely 3.3 km s"'. This value is, however, still well below the 
range of rotational velocities (vsin/ >8kms"') among which 
the frequency of Li-rich giants seems to increase. 

In RGB stars, the Li-rich phase has also been suggested to 
be associated with a mass loss-episode ( [de La Reza et al.[[1996| 
|1997| . This latter can cause blue-shifted asymmetric Hq, profiles 
or the presence of additional blue-shifted absorption lines in the 
Na doublet (NaD) U nes (Balachandran et al. 2000; Drak e et al.| 
2002[|Reddy & Lambert, ,200 5) in the stellar spectra. Figs. |6| and 



7J present a sample of the spectra obtained for the target stars 
in the Ha and NaD lines regions. Our stars do not present any 
evident asymmetry in the Hq, line (Fig.[6]l. #142173 is the only 
star that presents additional blue shifted lines around the sodium 
doublet (Fig.|7]i. These additional components are blue-shifted 
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Fig. 7. Sample of the target stars' spectra in the region of the 
NaD lines. 



by 40 and 39 km s^'with respect to the companion photospheric 
NaD lines. These velocities closely correspond to the local stan- 
dard of rest (vheiio=+38km s"'). Therefore, these features may 
originate in the interstellar medium, and not be of stellar origin 

[2000l l. 



(see, e.g., Balachandran et al. 



7. Conclusions 

We have built up a large database of high-resolution, high- 
quality stellar spectra of candidate thick-disk giants. We 
searched it for stars presenting a strong Li 1 resonance line and 
detected five Li-rich giants. One of them, namely #313132, has 
A(Li)>3.3 and therefore falls in the category of the rare super 
Li-rich giants (see Kraft et al. 1999 [Dom mguez et al.| 2004 
Monaco & Bonifacio , 2008 , [Kumar & Reddy 2009^|Carlberget 



al.| |2010| l 

All our stars but #313132 have kinematics compatible with 
thick-disk membership, and a comparison with theoretical evo- 
lutionary tracks (see Fig.|5]) of a proper metallicity is consistent 
with them being low-mass stars. They all seem to have evolved 
past the RGB-bump (Fig. [4]). Therefore, they do not belong to 
either of the two evolutionary regions where Li-rich giants tend 
to cluster, namely the RGB-bump or the early AGB in low- and 
intermediate-mass stars, respectively. Fig. |4] would suggest that 
these stars might have akeady undergone the extra-mixing pro- 
cess known to happen at the RGB-bump. The low carbon iso- 
topic ratios measured in #225245 may support this possibility, at 
least for this star. Under this hypothesis, they can be either pro- 
duced by CBP or are in the process of diluting Li again after the 
enrichment phase. In the latter case, their '^C/'^C isotopic ra- 
tios could be used as a diagnostic for the time scale involved in 
the dilution process after the Li-enrichment phase has occurred. 
#225245 has a quite high lithium abundance (A(Li)>2.7) and lies 
close to the RGB-tip. Its low '^C/'^C isotopic ratio is consistent 
with the CBP predictions of |Boothroyd & Sackmann| ( [T9"99j ). 



#313132 is a likely thin-disk member. Its spectroscopically 
derived atmospheric parameters are compatible with a ~4Mq 
star either in the RGB or AGB phase. It is the first super Li-rich 
giant detected at this phase (CBOO, GPBU09). Its '^c/i^c ra- 
tio is compatible with the first dredge-up predictions for similar 
masses. 

The detected Li-rich stars are slow rotators and do not show 
indication of mass loss, as shown by the inspection of the Hq, and 
the Na D lines. 
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